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^ - Abstract. An arbitrary initial magnetic field in an A star evolves into a stable 

equilibrium. Simulations are presented of the formation of non-axisymmetric 
equilibria consisting of twisted flux tubes meandering under the surface of 
the star, and analytic arguments are given relating the stability and form of 
these equilibria. These results may help to explain observations of Ap stars 
(~| , with very non-dipolar fields. This work is also applicable to other essentially 

non-convective stars such as white dwarfs and neutron stars. 
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. 1. Introduction 

Convective stars of various types tend to have small-scale, time-varying mag- 
^ ■ netic fields, which is interpreted as dynamo action driven by differential ro- 

^4 I tation and convection/buoyancy instabilities. In contrast, intermediate-mass 

^^O ■ main-sequence stars (1.5Mq < M < 6Mq), which have a small convective core 

and a radiative envelope, have either no detectable field or a large-scale, steady 
magnetic fields - the Ap stars (see e.g. Auriere et al. 2008). Many of these stars 
have roughly dipolar fields but many have rather more complex fields, such as 
T Sco (Donati et al. 2006). 
[ Given the lack of necessary dynamo ingredients, we infer that the magnetic 

fields are 'fossil' remnants in some stable equilibrium, formed from the mag- 
netic field left over from the convective/accretive protostellar phase. On the 
theoretical side, the emphasis has been on finding stable magnetic equilibria. 



I Analytically one can produce an equilibrium configuration simply by making 

sure all forces are balanced, and then test its stability using an energy method 
(Bernstein et al. 1958). Unfortunately, it has been easier to demonstrate the 
instability of various equilibria than to find stable configurations. For instance, 
both axisymmetric poloidal fields and axisymmetric toroidal fields have been 
shown to be unstable (Wright 1973; Tayler 1973; Braithwaite 2006). More re- 
cently, it has become possible to find stable equilibria using numerical methods. 
Braithwaite & Nordlund (2006, see also Braithwaite & Spruit 2004) found such a 
configuration by evolving an arbitrary initial magnetic field in time and watching 
it find its way into a stable equilibrium. This equilibrium is roughly axisymmet- 
ric, consisting of both toroidal and poloidal components (logical, since both are 
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unstable on their own) in a twisted-torus configuration. From the outside, the 
field looks approximately dipolar, as do the fields of many Ap stars. Here, I 
demonstrate the existence of more complex, non-dipolar configurations. I first 
describe simulations of magnetic equilibria evolving from arbitrary initial con- 
ditions, then use simple analytic methods to explore the nature and stability of 
these equilibria. 

2. Numerical simulations 

The setup of the simulations is described in detail in Braithwaite & Nordlund 
(2006); I give a brief summary here. The code used is the STAGGER CODE (Gudik- 
sen & Nordlund 2005), a high-order finite-difference Cartesian MHD code. The 
star is modelled as a ball of self-gravitating ideal gas of radius R arranged in a 
polytrope of index n — 3 (approximating an A star), surrounded by an atmo- 
sphere of very low electrical conductivity. The initial magnetic field resembles 
that expected at the end of a period of convection, i.e. chaotic and small-scale, 
containing wavenumbers with a flat power spectrum up to fcmax (set to 36R~^ 
here) , tapered so that B ^ fp where p is the gas density and p is a free param- 
eter. Note that if the star forms out of a uniformly magnetised cloud and the 
same fraction of flux is lost from all fluid elements, then we expect p — 2/3. 

The evolution in time of the magnetic field is followed in simulations with 
various values of the tapering parameter p: 0, 1/3, 2/3 and 1. In all cases, a stable 
equilibrium is reached after a few Alfvcn crossing times, but both the geometry 
and energy of the equilibrium is found to vary. Two types of equilibrium are 
found: simple axisymmctric (as in Braithwaite & Nordlund 2006) in the p = 1 
and p = 2/3 cases, and more complex non-axisymmetric in the p = 1/3 and 
p = cases. 

Fig. 1 shows the initial conditions (top- left), as well as both the p = 1 (top- 
right) andp — (bottom- left & -right) runs after around 9 Alfven crossing times. 
In the p = 1 case, the configuration is roughly axisymmetric and contains both 
toroidal and poloidal components. When p = 0, a non-axisymmetric equilibrium 
is reached. The configuration consists of twisted flux tubes below the surface of 
the star. The flux tubes lie horizontally and meander around the star at some 
distance below its surface. 

3. Stability of non-axisymmetric fields 

Consider a twisted flux tube of length s and width (at the surface of the star) 
aR whose axis lies at some constant radius in a star of radius R. The cross- 
section of such a flux tube is illustrated in Figure 2, as well as its axisymmetric 
relative. It is assumed that a <C 1. The three components of the magnetic field 
are toroidal (parallel to the length of the flux tube), latitudinal and radial; 
average values are denoted respectively by St, Bi and Bj.. The latitudinal and 



Non-dipolar magnetic fields in Ap stars 



3 




Figure 1. Top-left, the initial conditions. Semi-opaque shading represents the star; 
note that there is a lot of free space around the star, which ensures an accurate treat- 
ment of the potential (current-free) field outside the star. Top-right, the equilibrium 
reached in the case of a centrally-concentrated (p = 1) initial field. Bottom left and 
right (viewed from either side of the star) is the equilibrium reached in the p = case. 



radial components are collectively referred to as the poloidal component. Since 
the flux tube cannot have any ends, it must be joined in a loop, which can be 
wrapped around the star in any fashion. 

It is shown in Fig. 2 that the toroidal component of the field is confined 
to the volume enclosed by the largest poloidal line closed within the star, and 
is zero outside this volume. This has to be the case, because a poloidal line 
which crosses the stellar surface would 'unwind' if it contained any toroidal 
component. In this process, toroidal field is redistributed along the whole of the 
poloidal line. Toroidal field in the atmosphere is then destroyed by reconnection 
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Toroidal field confined 
to closed poloidal lines 



Axisymmetric equilibrium 




Figure 2. Cross-section of a flux tube in a star. In the axisymmetric case, the tube is 
a circle around the equator; in the non-axisymmetric case, the tube meanders around 
the star in a more complex fashion. Directions of Bi and Br are indicated by arrows 
towards the left; Bt is directed into the page. Poloidal field lines are marked with 
arrows; toroidal field is represented by shading. 

on a short time-scale, and the toroidal field inside the star is redistributed along 
the poloidal line to where it has been lost. This occurs on the Alfven time-scale. 
In this way, all toroidal field on the poloidal line is quickly transferred into and 
destroyed in the atmosphere. 

The magnetic energy of the flux tube is 



where s is the length of the tube (on the surface of the star), a is the angle 
subtended at the centre of the star. The three terms on the right hand side 
represent the field energy in the three components. The field energy outside the 
star can be ignored (Braithwaite 2008, in prep.). 

Allowing the length and width of the flux tube to change adiabatically, the 
tube will stretch or contract until dE/da — 0, i.e. an equilibrium is reached^. 
Since we are considering a magnetic field in a stably-stratified star where the 
magnetic energy density is very much less than the thermal, i.e. B'^ /Sir ^ P, we 
have the following restrictions on the displacements ^ which take place during 
this adjustment: 



^It is straightforward to verify that d^E/da^ > at this point, i.e. that the equilibrium is 
stable to stretching/contracting perturbations. 




(1) 



V • C « and g • C « 0. 



(2) 
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The consequence of these restrictions is that the volume of the flux tube stays 
constant since no matter can leave or enter it, which can be expressed by stating 
that the area of the flux tube on the stellar surface F = saR = const. In 
addition, we have: 

oina oma olna 

which follow simply from flux-freezing. Differentiating (1) with respect to a and 
setting dE/da = 0, and noting from the geometry that Bi w aBr, we have 

St « Bi « aBr. (4) 

Therefore the adjustment to equilibrium consists in stretching or contracting 
until the toroidal and latitudinal components are roughly equal. 

The toroidal and poloidal fluxes in the tube arc $t ~ (l/2)Q:i?,^J5t and 
$p PS {l/2)saRBr = {l/2)FBr respectively. They are conserved during adia- 
batic changes, and using (4) we have 

2 2 $t F $t ... 

so that if we know the area F and the ratio of toroidal to poloidal fluxes of a 
tube, we can calculate its equilibrium width aR and length s. 

The above arguments can be applied not only to the flux tube as a whole, 
but also to any small length of the tube with area 5F. Along the length of 
the flux tube <I>t is constant, but Br « 2S^p/SF can vary. A single flux tube 
need not have uniform width - at every point along its length, there is a local 
equilibrium giving a in terms of the toroidal flux and the local value of Br 
(remember from (3) that B,. is constant during adiabatic changes in a). 

In a star with more than one flux tube, neighbouring flux tubes must have 
their poloidal components in opposite senses to avoid boundary discontinuities. 
However, there is no constraint on the sense of the toroidal component, since 
the tubes are separated from each other by a region with zero toroidal compo- 
nent. Since magnetic helicity is essentially just the product of the poloidal and 
toroidal components, different tubes can have helicities of opposite sign, and 
configurations of zero net helicity are possible in principle. This contradicts the 
hypothesis that an arbitrary field should evolve into the lowest energy state for 
its helicity, rather the field can find local energy minima at higher energies. 



4. Discussion and conclusions 



Once an equilibrium has formed, it continues to evolve as a result of magnetic 
and thermal diffusivities {ij and k respectively). Magnetic diffusion (a result 
of finite conductivity) causes the field to evolve on a timescale C^/rj where C 
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is the relevant length scale. For an A star, this is around lO^^yr if £ = i?. 
Thermal diffusion can also cause field evolution. To be in pressure and density 
equilibrium with its non-magnetised surroundings, a flux tube must be at a 
lower temperature than its surroundings. Heat diffuses into the tube, causing it 
to rise on a timescale of order C^fi/n, where (5 is the ratio of thermal to magnetic 
pressure. These effects are discussed in Braithwaite 2008 (in prep.); it is likely 
that both are too slow to be important in A stars, except perhaps in stars with 
very complex, small-scale equilibria (i.e. with smaller C). Whether there is any 
observational evidence for field evolution in A stars is controversial. 

That we observe both axisymmetric and non-axisymmetric fields in Ap stars 
implies that initial conditions vary, and that probably the central concentration 
of the initial field varies. The reason for this is unclear. However, it is a common 
phenomenon that Ap stars share with white dwarfs and neutron stars (as is the 
enormous range in field strengths), while the apparent cutoff in field strengths 
below 300G (Auriere et al. 2008) is unique to A stars, so the two problems are 
probably unconnected. 

In summary, I have performed numerical simulations of the formation of 
stable magnetic equilibria in a non-convective star, starting from turbulent ini- 
tial conditions. A sufficiently centrally-concentrated initial field evolves into an 
approximately axisymmetric equilibrium; a more spread-out initial field evolves 
into a more complex, non-axisymmetric equilibrium consisting of twisted flux 
tubes meandering under the surface of the star. In this case, the adjustment to 
equilibrium consists in the stretching or contracting of these flux tubes until an 
energy minimum is reached. Still needed is a more quantitative understanding 
of the effect of initial conditions on the resulting equilibrium, and more work 
is also required on the quasi-static diffusive evolution of these equilibria. The 
results are applicable to all predominantly non-convective stars: main sequence 
above I.SMq, white dwarfs and neutron stars. 
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